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TECHNICAL MEMORANDUM 


THERMAL CONTROL OF HIGH ENERGY NUCLEAR 
WASTE, SPACE OPTION 

INTRODUCTION 


The space option for nuclear waste disposal presents some thermal 
problems which must be solved and assessed before feasibility by space disposal 
can be established. Primarily, the thermal problem is one of maintaining the 
waste temperature below the melt conditions. This is especially true for those 
waste compounds having very high energy densities ( >0.2 W/g) . Studies thus 
Far have resulted in development of the spherical waste configuration. Thus, 
the thermal problem is associated with the inner core of the waste package. 

The core temperature is 1 


T 


T + 
o 



4k 


( 1 ) 


The core temperature T is, therefore, greatly affected by the surface tem- 
perature, T^ . On the ground w here convection cooling can be made available, 

the core temperature is easier to maintain at acceptable levels. However, in 

space, the surface temperature is determined by the radiation properties of the 

waste surface. The surface temperature is established by its area and 

emittancc. Increasing the surface area decreases the surface temperature but 

increases the second term, q R 2 /4k. The net result is always an increase in 

o 

the core temperature with an increase in R . 

o 

Thus, in space, very little control can be exercised over the core 
temperature. For a homogeneous waste distribution, the temperature is best 
limited by limiting the energy generating characteristics and mass size. How- 
ever, this is self-defeating since nuclear waste has an inherently low conductivity 
and some of the most desirable waste to be disposed of has high energy content. 
Also, it is desirable to use the Shuttle payload capability for each launch. 

J. Kreith, Frank: Principles of Heat Transfer. Dun-Donnelley, New York, 1976. 



The purpose of this study is to suggest a technique whereby the effec- 
tive conductivity of the waste can be increased. This allows core temperature 
control with larger masses having larger energy contents. The concept is to 
bring the energy within the core to the surface by means of cylindrical rods 
having high conductivity. This does not reduce the surface temperature but 
does reduce the temperature difference occurring between the surface and core. 

The objective herein is to evaluate this concept and to establish the 
effectiveness of the cooling rods as greater waste masses having greater energy 
densities are considered. Also, the effects of the number of rods and their 
size will be assessed. 

The cooling rod concept is also important because of the way it lends 
itself to manufacture of the nuclear waste package. This advantage is very 
important and is the motivation for studying the rod concept in detail. This 
cooling rod packaging process consists of three basic steps (in order) : 

1) A steel hemispherical shell with rod holes located on its surface, 
which will later receive the rods, is packed with waste material by a suitable 
packing mechanism. There are thin membranes over the rod holes to contain 
the material during the packing process. 


i PACKING 
I MACHINE 
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2) After packing two hemispheres, they are brought together and 
welded. The cooling rods are then inserted through the rod holes. There may 
be as many as 300 to 500 rod holes. The cooling rod pierces the membrane as 
it is driven into the waste material. 



3) After each rod is driven flush, it is welded to the sphere's 
surface. The package is now complete for deportation into space. 



There are several advantages to this technique besides the symmetrical 
configuration which lessens the thermal analysis problem. First, there is less 
chance of waste spillage and tool contamination. This is an important advantage 
and motivates use of this concept. Second, the extra packing of the waste 
material as a result of driving the rods reduces the contact resistance between 
the rod and waste material. 
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GENERAL THERMAL CHARACTERISTICS 


The temperature distribution within a spherical waste package having 
a homogeneous mass distribution can be found in Footnote 1* 



The temperature, T, occurs at Radius R, with R q being the radius of the 

spherical package. At the core, R/R q equals zero and equation (2) reduces to 

equation (l). At the surface, R/R q is unity, and the temperature in question 

if equal to the surface temperature is determined by the radiation characteristics 
at the surface: 



Equation (3) states that in the steady state the total energy rate generated must 
be radiated at the surface. The required surface temperature is T q . 

Equations (2) and (3) are very simple to apply but are not in terms of 
the variables which the nuclear material is usually specified. In general a waste 
mass, M, is given which has a given energy density, E, and a mass density, p. 
Noting that: 


q = Ep 


R = 
o 



v s 


Q = EM . 
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Equation (2) can be manipulated to give a more useful form in terms of stated 
characteristics. These substitutions in equation (2) give 



A typical temperature distribution resulting from this equation is given in 
Figure 1. In this particular example the energy density of 0.002 W/g is 
sufficiently low that cooling conduction rods are not required. 

Most of the time, only the core temperature is of interest since 
therein lies the highest temperature. Additional characteristics for just the 
core temperature are given in Figure 2. Waste mass is the independent vari- 
able with energy density as an argument. A density of 2. 88 g/cm 3 and con- 
ductivity of 1 s \V m-" K arc taken as typical waste characteristics. It is noted 
that density and conductivity are illusive parameters. The conductivity is 
particularly difficult to characterize. 

Figure 2 illustrates the difficulty of keeping the core temperature 
below 1800°F for energy densities above 0.004 W/g. Yet, some nuclear waste 
lias densities above 0.2 W/g. Cooling rods are necessary for these energy 
levels. 


Figure 3 has been included to illustrate the importance of conductivity. 
This figure emphasizes the importance ot achieving an effective conductivity of 
at least 4 or 5 W/m-°K. From equation (4) , it is realized that as the conduc- 
tivity becomes greater, the core temperature approaches the surface tempera- 
ture. The objective herein is to demonstrate how the cooling rod increases the 
effective conductivity of the waste material. 

The effective conductance can be found bv solving equation (4) for k 
and taking the value of T as the core temperature which results from application 
of the cooling rods: 



TEMPERATURE 9 R — °F 



Figure 1. Temperature distribution for a spherical waste configuration having 
a homogeneous mass distribution. 




CORE TEMPERATURE - °F 


DENSITY 2.88 g/cm 3 
CONDUCTIVITY 1 .8 wlm • °K 



2 3 4 5 6 7 8 


WASTE MASS - kg X 10 3 


Figure 2. Core temperature as computed from equation (4). 




( 5 ) 


Equation (5) will be utilized later to compute the effective waste conductivity 
resulting from cooling rods. 




5 1.0 15 25 

CONDUCTIVITY - (W/m °Kt 

Figure 3. Core temperature sensitivity to nuclear 
waste conductivity. 


Equation (4) has been computerized for quick assessment of core 
temperature. Also, if the effective conductivity is known, the core temperature 
can be quickly computed without resorting to the SINPA system numerical 
difference analyzer program. A typical output from this program is given in 
Table 1. The program iterates waste mass up to 10 000 kg and energy densities 
up to 0.008 W/g. The values of these selected ranges can be easily changed as 
is discussed in Appendix A. In some material, the temperature which can be 
tolerated is not known. For this reason, a tool for quick evaluation of any pro- 
posed nuclear waste package can be made available to those whose expertise are 
with nuclear material. 
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THERMAL MODa GEOMETRY 


As already mentioned, the core temperature cooling mechanism 

consists of conductive rods which penetrate the waste radically within a few 

Inches of die spherical center. The rod radius, r , results in interference 

o 

between rods as they approach the center (Fig. 4) . The radius of the resulting 
sphere at the center can be found by equating the area of the inner sphere with 
the sum of the cross-sectional areas of the rods plus 10.7 percent: 

4*r 2 = 1.107 rNr 1 
o c 

r_ 

r = -r 2 V 1.107 N 
c 2 


For 362 rods, the inner core diameter will be 2. 5 in. for 0.25-in. radius rods. 
This core will probably be made of steel to provide structural integrity to the 
waste package, in any event, the allowable rod length must be accounted for by 
this method. The first step to be taken is to define a relationship for the number 
of rods. At first, it may appear that any arbitrary number could be used. 
However, this is rot the case if the rods are to be equally distributed through- 
out the sphere. 

It is necessary to establish the relationship between the number of 
rods and the distance between rod tips as they enfr*r the spherical surface. 
Finally, the distance between all adjacent rod tips would be equal. The basis 
for establishing these relationships is a surface constructed of 20 equilateral 
triangles. The solid figure formed by a single triangle is known as an 
icosahedron. Representation of a spherical surface by 20 icosahedrons is 
shown in Figure 5. 

Each equilateral triangle can be- subdivided into small triangles as 
i' .ustrated in Figures 6 and 7. Each intersection lies on the surface of the 
sphere. As the number of triangles increases each icosahedron surface 
approaches an equilateral spherical surface. In Figure 6 the basic icosahedron 
leg is divided equally into two elements. This results in formation of four 
triangles. The number of equal subdivisions in each icosahedron sides is known 
as the geodetic frequency. F igure 7 has a geodetic frequency of six. 
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r igure 4, Cooling rod concept which conducts heat energy from the 
core to radiating surface ol the waste. 

For purposes of the thermal control mechanism, a rod would be 
located at each interaction. The relationship between the geodetic frequency 
and the surface characteristics arc: 

Number of Vertices 10 r 2 + 2 

Number of Faces 20 v 2 


For example, a geodetic frequency <4 h results in 3h2 rods and 720 triangles 
formed by each rod end. Since the frequency must be an integer, the possible 
number of rods grows as 12, 12, 02, 102, 2f>2, 302, 402, etc* Ail intermediate 
numbers of rods art not possible with this configuration. 
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CO 


Figure 6, A spherical surface generated by a geodetic frequency of 2, 





Figure 7. A sphere al surface generated by a geodetic frequency of 6. 


For a frequency of one* the geometry reduces to the basic 20 equi- 
lateral triangles. Sides which form the surface of the icosahedron are equal. 
However* geodetic frequencies above one do not produce all equilateral 
triangles. The subdivided triangles do not result in equilateral triangles, 
although some sides are repeated. Consider the icosahedron surface of Figure 
8 where the sides have been subdivided with a frequency of four. Letters indi- 
cate members of the same length. The length of each member can be found by 
multiplying the encompassing sphere radius by a correction factor* CF. Some 
of these correction factors are given in Table 2. For an example, the length 
of the "c" member for a frequency of 4 is the radius times 0. 294530. It is 
unfortunate that this situation exists; however* for purposes of the thermal 
model and rod placement, the problem becomes greatly simplified by generaliz- 
ing on a single relation between the sphere radius and element length. For 
purposes of the thermal model, the following generalization has been made: 


R 

o 


l 


0. 8684 
v 


( 6 ) 


Equation (G) is based upon the length trends occurring within the icosahedron. 

It is noted that this equation cannot be used for actual design construction but 
docs allow thermal analysis of a single rod to represent all rods. 

Figure 9(a) is a plain view projection of a portion of an equilateral 
spherical triangle. All of the element lengths are equal as per equation (6). 
Each intersection is the location of a rod. Around each rod is a prescriber 
circle of diameter, f , the clement length. Each of these circles represent a 
cone, each having identical thermal distribution. 

An element cone is shown in Figure 9(b). This cone is defined by 
dimensions ( and R^ as related by equation (G) . The number of cones is 

equal to the number of vertices which is, 10 r 2 + 2. Thus, the mass element to 
be thermally modeled is defined. The complex thermal and mechanical config- 
uration has been reduced to a simple cone w ith a single rod located on its radius 
centerline. The thermal distribution in this cone characterizes the thermal 
characteristics of the entire nuclear waste package. It is noted that a small 
area exists outside of the prescribed circles of Figure 9(a). 

The energy represented by this area is distributed proportionally 
w ith in each cone. Thus, the temperature at the outer surface of the cone is 
considered to be representative of the temperature between intersecting cones. 
This deviation from the actual configuration is considered to be minor and the 
simplicity it brings to the model is warranted. 
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Figure 8. Subdivision of an icosahedron to form an equilateral spherical surface. 
TABLE 2. ELEMENTS LENGTHS, ( , CORRECTION FACTORS (CF) 































THERMAL MODEL NODE LOCATIONS 


Now that the mass element to be thermally modeled has been defined, 
the node location and their respective capacitance and conductance are to be 
determined. Since each node repress nting the nuclear waste is an energy 
generating node, it became an intuitive guideline to give each node equal mass. 
In this way, the energy generated would be distributed equally within the model 
network. This approach reinforces the fidelity of the model since a finite 
number of nodes represents an infinite number of energy generating points. 

To accomplish this guideline the sphere was divided into four concentric shells 
each having the same mass. This subdivision of the element cone is illustrated 
in Figure 10. Each shell is numbered 1 to 4 and they cut through the cone as 
shown. For proper location of the shells, each section of the cone will have the 
same mass. If n is taken as the shell number, the necessary relationship 
between the first shell, n= 1, radius, R , and the radius R^ , at n is 





3 

1 


R 

n 



(V 


The spherical radius, R , will be taken as the independent variable with n = 4. 
Therefore, R^ becomes 


R 


2—^ — = R 
1/3 1 


( 4 ) 


where 


( 8 ) 


R 


n=4 


R 

o 
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Figure 10. Spherical shells having equal mass. 

For an example, if the waste radius is 2.4 ft, R^ becomes 1.5119 ft. From this 
value equation (7) can be used to find R^ at other shell locations. A scaled 
sketch of the waste package divided into shells as above is given in Figure 11. 

To complete the thermal model, the mass within the top three shells 
is divided into doughnuts as shown in Figure 12. The generalized radius of the 
outside doughnut is r ^ . The inside doughnut radius is r . For equal mass 

within each doughnut, the following relationship exists: 




( 9 ) 


For similar triangles. 


1/2 

R 

o 


r 


n2 


R 

n 
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Thus 


R 

r n2 = 2(R ~7t) 


which determines r of equation ( 9) . Equations (7) through ‘ ’0) completely 

describe the geometry of the mass element. A summary of these dimensions 
are given in Figure 13. 



Figure 13. General dimensions of the mass element. 
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The nodal system consists of six doughnuts, each having different 
dimensions as given in Figure 13. Also, there is node 10 which represents a 
cone shape. The mass in the cone is twice the mass in any single doughn ut. 

The cooling rod common to each set of doughnuts and cone has nodes 50, 40, 

30, and 20, respectively. The R-C network representing the model is illustrated 
in Figure 12. The nodal system consists of 13 nodes with 99 being a space node 
and node 17 the plate which contains the waste material. 

It is noted that the model as described here has not been presented to 
tiie depth where as equations are available for computing the capacitance and 
conductance between nodes. The intention is to describe the model in sufficient 
detail that the modeling approach can be understood. Since it was anticipated 
that the model would have to be integrated many times, the model was completely 
computerized and the output was made compatible with the SINDA data block 
forms. A typical output from this program is given in Figure 14. This pro- 
gram is discussed in Appendix B. 


DATA RESULTS 


Before the computerized thermal model was complete, it was diffi- 
cult to speculate upon the effectiveness of the cooling rods. The entire purpose 
was to assess the capability of the rods in controlling the core temperature. 
Results from one of the first sample cases are shown in Figure 15. While the 
specific temperatures obtained from this sample problem are of an academic 
interest only, it indicatec that the cooling rod concept can be very effective in 
reducing core temperature. For this sample case, a mass of 10 480 lb was 
chosen with 362 half inch rods. The energy density was 0.2071 W/g. Without 
the cooling rods the core temperature was ove 7 ' 80 000°F. With the cooling rods 
the core temperature dropped to 8500°F. This is a factor of 10, a result worth 
noting. The surface temperature is about 1900° F. The surface temperature is 
the same with and without the rods since, in the steady state, the same energy 
must be radiated at the surface. The cooling rods used in this example were 
molybdenum which has a conductivity of 72 Btu/hr-ft- # F at 3000 °F. 

To establish the limit of capability of the cooling rods, the passive 
molybdenum rods were replaced with heat pipes. These results are indicated 
in Figure 16. Most of these temperatures can be tolerated. Thus, in the limit, 
the rou cooling concept can maintain core temperature for large masses and 
high energy densities. It is of interest to note that the reduction in core tem- 
perature with the number of rods is not significant. That is, the reduction in 
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179. MM 


UASTE DENSITY LB/TT3I 
RADIUS OF CONDUCTING ROD l PT 
GEODETIC FRE9UENCV 
HASTE RADIUS! FT 
* HEAT CAPACITY, WASTE i »TU/X» F 
HEAT CAPACITY, ROD l DTU/LD F 
CONDUCTIVITY, WASTE! DTU/I9I-FT F 
COHDUCTIWITV, ROD! DTU/HR-PT F 
ROD DENSITY! LB/FT3 
E NERGY DENSITY! BTU/MI LB 
NUMBER OP VERTICES 
TOTAL NASS! LBN 
RADIUS OP FIRST SHELL! FT 
RATIO OP GEODETIC ELEMENTS 



1B4M!b3U 


1.S11S 

6.SS73 


BCD 3N00E DATA 

11. 2M., .4459 

12.2M.. .4459 

13.2M., .4459 

14.2M., .4459 

1S.2M., .4459 

16, 2M., .4459 

17,2M., .9874 

19.2M., .8917 

2B.2M., .9536 

39.2M., .92M 

49,2M., .9149 

S9.2M., .9111 

DCD 3SOU0CE DATA 

11, 648.5378 

12, 648.5378 

13, 648.5378 

14, 648.5378 

15, 648.5378 

16, 648.5378 
19,1297.9755 

DCD 3C0HDUCT0R DATA 

1112.11.12. .3913 

1213.12.13. .7924 

1415.14.15, .4995 

1516.15.16, .7217 
1114,11,14.12.3856 

1215.12.15, 8.6442 

1316.13.16, 6.8675 
1911,19,11, 1.9529 

1914.19.14. 1.1136 
1929,19,29,11.5998 

1439.14.39, 4.31M 

1549.15.49, 3.9295 

1659. 16. 59, 2.4118 

1317.13.17, 1.9244 

1617.16.17, 1.9244 

5917.59. 17, .5729 

2939.29.39, .9659 

3349.39.49, .1877 

4959. 49.59. .2535 
-1799.17.99, 1.673189-19 


Figure 14. Typical computer output which generates 
nodal data for SINDA. 
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core temperature by using 810 rods is not significant from that of 492 rods. 
This is especially true for waste masses below 10 480 lb. Thus, geodetic 
frequencies greater than 7 (492 rods) is not warranted in view of the additional 
complexity of building the waste package. The effective conductivity of the 
entire package was calculated for the 10 4S0-lb package u. ng equation (5) . 

The effective conductivity was computed at 122 W/ra-'K. 

It became of particular interest to assess the performance of the 
cooling rods by reducing the waste mass while keeping the same high energy 
density. These results are shown in F igure 17. Three cases are shown for 
different diameters of rods. Since most waste temperatures must be less than 
3500*F, the waste mass must be limited to 2000 to 3000 lb. This implies that 
waste material which has high energy density must consist of multiple packages 
of relatively low mass. 

Figure 18 corresponds to Figure 2 except that much higher energy 
densities are considered as the argument. To approach a single waste package 
of 10 000 lb, the energy density must be below 0. 1 W /g to keep the core tem- 
perature below 3500°F. 
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Figure 18. Core temperature with energy density as an argument, 


CORE TEMPERATURE - °P 
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APPENDIX A 


TEMPERATURE DISTRIBUTION FOR A SPHERICAL 
HOMOGENEOUS MASS 


The computer program presented in this appendix represents the 
solution to Equation (4). There are three inputs: waste density, waste con- 
ductivity, and surface emittance. These parameters are noted as RHO, XK, 
and EMISS, respectively per statement 17 of the main program. Density has 
units in g/cm 5 and conductivity in W/m-°K. An emittance value of 0. 80 has 
been selected for all calculations herein. The range of mass is controlled by 
state 19 of subroutine PRNT. As shown the mass is iterated in increments of 
2000 kg, beginning with 1000 kg up to 10 000 kg. 

For each mass the energy density is initiated from 0.001 to 0.005 
W/g as controlled by statements 44 through 48 of the same subroutine. The 
mass and energy values can be alternated to achieve any combination of mass 
and energy density desired. It is noted that SIG, the Stefan Boltzmann constant, 
has units of W/m 2 -°R. 

Vol — Volume - m 5 

RADI — Radius — m 

QDOT — Energy Rate — W/m 3 

QT — Total Energy Rate — W/kg 
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S4* 

Site 
SSI 
S 3 1 

sorts? scuNiit 


t#3 forna? a#K,3SH surface enittiince # Fi§*i # ///t 

m ?wmm c t2K # 4HC0S€ f i4K,4NeoRc, m* 

its ?mmt t tm,wTtm, lm L m mm, im> mmm, i3K, tmm* m» smmm* 

tm rmm n t i 

ilEK, 3HCF^ f I IK, TWMRTTS X liK ( 9MCy*?T$43, /✓ 1 
m format .ri?a,riia # ri?a # il?,ri#,3 # wa,rii.4? 

RETURN 
END 


FCRROR 27R3AN4 C33 $l?3Rt §4/|?/?9 MMMMI4 

jwt?>Mnam*miu%t %r.mm 

subroutine c mu 

confutes equations for mmm 


i 

a 

3 

4 

5 

6 
7 
§ 
3 

it 

II 

ia 

13 

14 

15 
IS 
I? 
IS 
13 
20 
21 
22 

23 

24 
26 
28 
27 
21 


C 

c 

c 

c 

c 

c 

c 

c 

c 

€ 

€ 

c 

c 

c 


IHCUHME CORK* 

wwimus mz m roiious 
mt ~ mum 
nmt ~ mmm 

mm ~ mum mm mu mum 

a? - tmm mm 

sis * swfm-minmm constant 

TO - SURFACE TEMPERATURE 

re - surface - mmn 

t - core mmMmvm ~ mwtn 

T f - CORE TERPERRTIRE - FAHRENHEIT 

U0L * C R/RNO III , §€>#3 

tA0I*C3. tWUt 4.13, 14 1 lit, 333 

0D0T*EfRHO/| ,§E~#S 

QT*i*WMESR 

$I0*5*407£-09 

T0*C0T/{4*t3* 1 4 tRRDl if2fSISi£RlS$? >S* * 26* 

tf* c to- 32, mw-m* mn. is 

T»TE4ai9TjmADIti 2/C 4 , f XJC 1 
TF*9.ItT-2?3*16 J/5.432# 


4S#* 


RETURN 

END 



APPENDIX B 


PRESINDA, SINDA, DATA BLOCK DEVELOPMENT 


This appendix presents the PRESINDA program. This program 
computes the conductance, compacitance, and energy rate of energy generating 
nodes. There are 10 inputs as noted in statement 24 of the main program. 
Units are given in SUBROUTINE PRNT statements 136 through 145. 


32 



JHEERS1IN2888JPPILE9U KPRES1NDA 
. 1 C ttPRESXNOA** 

2 C 

3 INCLUDE COflDCK 

4 C 

5 C INITIALIZE TEKTRONIX GRAPHICS PACKAGE 

6 CALL INITTU29) 

T CALL TERN) 3. 1924) 

8 CALL CMRSIZO* 

9 CALL TSEND 

19 C 

11 REWIND 8 

12 C 

13 C LOOP - UNTIL ALL CASES HAVE BEEN DONE 

14 C RHU - UASTE DENSITV 

1C C RO - RADIUS OP CONDUCTING ROD 

16 C XNL - GEODETIC FREQUENCY 

17 C R JASTE RADIUS 

15 C CPU - HEAT CAPACITY. JASTE 

19 c CAR - HEAT CAPACITY. ROC 

£9 C XKU * CONDUCTIVITY, UAS T E 

£1 C X<R - CONDUCTIVITY, ROD 

cB C RHO - ROD CENSITv 

£3 C ED - ENERGY DENSITY 

£4 19 READ :S, 199. END-999) RHU, RC. XNU.fi .CPJ.CPR, XKU. 

25 tXKR.RHO.ED 

26 C 

£7 C DO COPP JT AT IONS 

28 CALL COUP 

29 C 

39 C PRINT OJT DATA 

31 CALL PRNT 

32 C 

33 C LOOP - BACK TO NEXT CASE 

34 GO TO 19 

35 C 

36 199 PORPAT « ) 

37 999 STOP 

38 END 


oo 

w 



00 

.u 


11 

a*c 


SUBROUTINE PRNT 

3* 

4ic 


INCLUDE COMDEX 

SI 

6tC 


CALL NEUPAC 

7»C 


PRINT INPUT DATA 

8* 


PRINT :el,RHU 

91 


PRINT 1D2.R0 

18* 


PRINT 103. XHU 

1U 


PRINT 104. R 

12 t 


PRINT 1 05. CPU 

13 * 


PRINT 106. CPR 

14* 


PRINT 107, XKU 

15 * 


PRINT 108. XKR 

16* 


PRINT 109. RHO 

17* 


PRINT 110. ED 

18* 


PRINT lll.XN 

19* 


PRINT 112.XP 

26: 


PRINT 113. RIR 

21 * 
22 : C 


PRINT 114. RL 

2HC 


PRINT NODE DATA HEADINGS 

24* 


WRITE (8,115) 

25* 


WRITE <8,117) 

26* 


URITE (8,118) 

27t 


WRITE (8,119) 

28* 


PRINT 119 

29* 
36 * C 


FD«200. 

31 :C 


SET NO EQUAL TO NUN8ER OF NODE 

32 * C 


FOR USE IN UARIA0LC FORHAT 

33: 


ND»10 

34* 


DO 10 N-1.7 

35* 
36 *C 


ND«ND*1 

37 1 C 


WRITE NODE DATA TO A 

38 * C 


FILE AND ON TO SCREEN 

39* 


URITE (8,140) ND.FD.DN(N) 

46* 
41 * C 

16 

PRINT 140,ND,FD,DN(N) 

42 *C 


SET ND EOUAL TO NUMBER OF NODE 

43 *C 


FOR USE IN VARIARLC FORMAT 

441 


ND*0 

4St 


DO 20 N»B. 12 

46* 

4?*C 


ND-NDM0 

46* C 


URITE NODE DATA TO A 

49* C 


FILE AND ON TO SCREEN 

St> 


WRITE (8,140) ND.FD.DN(N) 
PRINT 140.ND.FD.ON(N> 

91* 

86 

58* 


URITE (8,141) 

S3: 
64 1 C 


URITE (8.118) 

SSIC 


PRINT SOURCE DATA HEADING 

66* 


PRINT 120 

671 


URITE (8,120) 


SS»C 

S9*C URITE SOURCE MTR TO A 

66*C FILE RHD OR TO SCREEN 

Sli DO 30 IS* 1 ,6 

Ut URXTE <8,t48> IS,StIS> 

63* 3# PRINT 142.IS.SCIS) 


64 * C 
6S*C 

SET IS EOUAL TO NUN8CR OF NODE 

66 *C 

FOR USE IN A VARIABLE FORMAT 

67 « 

ISMO 

68* 

URITE (8,143) SC7> 

69* 

PRINT 143, SC7) 

76* 

URITE (8.118) 

?1*C 
72 1 C 

PRINT CONDUCTOR DATA HEADING 

731 

URITE (8,121) 

741 

PRINT 181 

75 *C 

76 *C 

SET IC AND JC EOUAL TO 

77 »C 

NUMBER OF CONDUCTOR NODES 

78 *C 

FOR USE IN UAPIA8LC FORMAT 

79* 

ICt 1 )*11 

86* 

JC(1 >-12 

8X i 

IC(2 >*12 

82* 

JC (2 )• 13 

83* 

ICO >»14 

84 t 

JCO)-IS 

85 1 

IC(4)-1S 

86* 

JC ( 4 )»16 

87t 

IC(S)-11 

88* 

JCC5J-14 

89* 

IC (6 )■ 12 

96* 

JC (6 )• IS 

91* 

ICC?>»13 

92* 

JC(7)-16 

93* 

IC<8 )*10 

94* 

JC(8)-11 

95* 

IC(9)-10 

96* 

JC(9)*14 

971 

IC(10)*10 

98* 

JC( 10 )*20 

99* 

ICC11 )-14 

166* 

JCdt )«30 

161* 

IC(12)-1S 

162* 

JC< 12 )"40 

163* 

IC(13)*16 

164* 

JC< 13 )*S0 

165 * 

ICC 14 >*13 

166* 

JC( 14 >»17 

167* 

IC( 15 >*1S 

166* 

JC< 15 )«17 

169* 

IC(16)*S0 

116* 

JC(15)*17 

111* 

IC( 17 )»80 

118* 

JC(17»30 

113* 

IC( 18 )*30 

114* 

JC( 18 )«48 


U4i> 



1171 
1101 
1 19*C 
12B*C 
121 *C 

122 i 

123i 

1241 

1 25 1 

126* 

127t 

128*0 

129*0 

1 30i 

13UC 

1 32*C 

133* 

134: 



178* 

142 FORMAT 

CtlX 

IC< 19 > *4# 

173* 

143 FORMAT 

<ux 

JCt 19 >*50 

174* 

144 FORMAT 

U1X 

1C<20 )*17 

175* 

145 FORMAT 

tux 

JC< 20 >*99 

) 

176* 

146 FORMAT 

('88: 

WRITE CONDUCTOR DATA TO 

177*C 



A FILE AND ON TO SCREEN 

178* 

RETURN 


DO 40 >1,19 

179* 

END 


WRITE <8,144) ret I ), f C( n, ICt t >, JCtlj.CC I ) 

EOF: 179 SCAN *64 



40 PRINT l44,IC<n.JCClJ.XC<X>.JCcn.C(I) 0O 

WRITE <8,145) IC(20>, JC(20),IC<20),JC<2e>.f <20> 
PRINT 145, XCt 20), JCC20), IC<20>,JCt20 >,C<2* < 

WRITE (8,118) 

PRINT I'S TO SEPARATE DATA 
CASES WRITTEN ON FILE 
WRITE <8,145) 


1,'.',F9.4> 

,F«.4> 

,'*'.12.', '.12, 

s.xa.'.'.ia.',' 

8t8t ' ^ 


CALL WDCOPV 
CALL NEUPAG 
CALL TSEND 


1 35 ? C 

136: 101 FORNAX 
137t 102 FORMAT 
138* 183 FORMAT 
139 i 104 FORMAT 
1401 IBS FORMAT 
1411 105 FORMAT 
142 * 107 FORMAT 
143* 108 FORMAT 
1441 109 FORMAT 
145* 110 FORMAT 
146* 111 FORMAT 
147* 11C FORMAT 
148* 113 FORMAT 
149* 114 FORMAT 
150* 116 FORMAT 
151* 117 FORMAT 
152* 118 FORMAT 
153* 119 FORMAT 
1541 120 FORMAT 
155 * 121 FORMAT 
156*C 122 FORMAT 
1S71C 123 FORMAT 
158 * C 124 FORMAT 
159*0 125 FORMAT 
160 I C 126 FORMAT 
161*0 187 FORMAT 
162 <C 128 FORMAT 
163 iC 129 FORMAT 
164*0 130 FORMAT 
166*0 131 FORMAT 
166*0 132 FORMAT 
167*0 133 FORMAT 
161*C 134 FORMAT 
169*0 135 FORMAT 
176* 140 FORMAT 
1711 14X FORMAT 


U0X,35H WASTE DENSITY LB'FT3* 

<10X,3SH RADIUS OF CONDUCTING ROD* FT 
<10X,3SM GEODETIC FREQUENCY 
( 10X, 35M WASTE RADIUS* FT 
<10X,3SH HEAT CAPACITY, WASTE* BTU/Li F 
<16X,35H HEAT CAPACITY, ROD* BTU'LB F 
< 10X, 35N CONDUCTIVITY, WASTE * BTU'HR-FT 
<10X,3SN CONDUCTIVITY, ROD* BTU/HR-FT F 
<10X,3SH ROO DENSITY* L0/FT3 
<10X,3SH ENERGY DENSITY* 0TU/HR LB 
U0X,35N NUMBER OF UERTICES 
OBX,3SN TOTAL MASS* L0M 
aex,35H RADIUS OF FIRST SHELL ! FT 
U0X,33M RATIO OF GEODETIC ELEMENTS 
(6X, 18H BCD 3THERMAL SPCS ) 

1 6X, 12H BCD 9NEWMAN) 

<6X,4H END) 

(6X,15H BCD 3NCDE DATA) 

<6X # 17H BCD 3S0URCE DATA) 

(6X,20H BCD 3CONDUCTOR DATA) 

(6X,2#H NOD 3C0NSTANTS DATA) 

(9X.12H 1 1 MEND* 120 , > 

<9X,16H OUTPUT -6. ) 

< 6X, 16H BCD 3ARRAV DATA) 

(6X, 15H BCD 3EXECUTI0N ) 

( 22WF DIMENSION X<200)) 

U4MF NDXM*206> > 

< 1 1HF NTH*# ) 

(9X.7H CNFRWD ' 

<6X,17H BCD 3UARIABLES 1) 

(6X.17H BCD 3 VARIABLES 2) 

C6X.16H BCD 30UTPUT CALLS) 

(9X.7H TPRINT) 

C6X.17H BCD SEND OF DATA) 

uix.ia. ',F4.t,' # %F7.4> 
tux, 13H-99, -456, , . 06 ) 


, r 20.4> 

, F20 • 4 ) 

, r 20 . 4 ) 

, r 20 • 4 ) 

, FB0 . 4 ) 
,^20. 4) 

F , Fg0 . 4 ) 
, r 20. 4) 
,F20.4) 
,F£0.4) 
,^20.4) 

« F20 . 4 ) 
,r2e.4 ) 
,*20.4,0 


' » ' *F7. 4) 
,I2.','.lPEia.6 



*70 


w 


SCI MOOT INC CO 


INCLUDE CONDCK 


COMPUTE EQUATIONS POP VARIAOLCS NOT READ IN 
XN - NUNOER OP VERTICES 
XN - TOTAL NASS 


RIP - I .’AD I US OP PI® 
RL - RATIO OP CC . 


SHELL 

IC ELEMENTS 


XN»iA.«XNU»«a>2. 

XN* l 4 . 1 B 878 Ra 83 -C 3 .)-» 8 R 0882 aR 8 XN/RM 0 ) )*PHU 


RIR*R/4.88<l./3. > 


RL*XHU8a.B6B« 


SET VARIABLES EQUAL TO EQUATIONS USED IN COMPUTING 
NODE. SOURCE, AND CONDUCTOR DATA EQUATIONS 
Al*<2.8((2./3. ))8<RIfi882)/<4.8< RL882 ) ) 
A2*(2.88<2./3. ) >8<RIR882>/<8.8<RL8a2) > 
A3*(3.88(2./3. ) )8(PIR882 )/(4.8«RL8a8> > 

44* ( 3.88<2./3. ) >8<RIR888 )/(8.8<RL888 > ) 
AS*<RIRaa2)/<8.a<PL882>) 

AS* ( RIR882 )/< 4. 8< RL882 ) > 

A7*(<4.88(2./3. >>8(R1R888) >/<4.8(RL882>> 

A8*( (4.88C2./3. >>8(RIR882)>/(8.8<RL882n 
1 1 • ( ( 3 ■ ti ( l , /3 . >>8RIR> 

12* < < 4 . 88 ( 1 . /3. ) >8RIR ) 

B3* < (8. 88 C 1 . /3. >>8RIR) 

84*RIR/(4.8RL) 

BS-< <2. aac i ./3. i)8RiR>/<4.8RL) 

BS*< o.aaci./s. > >8RiR>/(4.aRi > 

B7*(Roaa2)/a. 

88*0.14/4. >a((2.aR0)aa2) 


COMPUTE EQUATIONS FOR NODE DATA 

SET NODES II THRU IS EQUAL TO SAME EQUATION 


DC 18 N* 1.6 

10 Oh<N>*xNa.?saCPU/<XNa6. > 8 

DN<7>*3.148A781?8.8.S28.2 8 

DN<8>*xna.25aCPU/XN 8 

DN(9 »*3.148<ROaa2>a<RIR-$ORTO.aXN>a 

a((2.aR0>/4. naRHoacPR 8 

DN(10>*3.14a( R0882 >8 < 83-RI R ) a RHO 8CPR 8 

DN( 1 1 )»3.l48(R0888>8(tl-83>8RM08CPB 8 

DN(i2>*3.i4a(Roaa2iai82-8i isRHoaCPR f 


COMPUTE EQUATIONS POP SOURCE DATA 

SET NODES 11 THRU IS EOUAL TO SAME EQUATION 

DO 20 IS* 1 , 6 

28 S( IS )*NMaEDa.7*/(XN*8. ) 8 

S(7)*XN8ED8.8S/XN 8 

COMPUTE EQUATIONS FOR CONDUCTOR NODES 
C (I ) *3 . 1 48 ( A1 -A2-87 )/ C . sa ( * 1 -RIP ) >8MCU Q 

C< 8 > -3 . 148 ( A3-A4-87 )/ ( . 88 ( 82*83 > )8MCU • 

C<3)*3.148tA2*BT)/( .Sa ( ll-RIR ) )8XKW Q 

C< 4 )*3. 148 < A44B? >/l .58(82-83) >8Xru Q 

C<S>*<2. 8(3,14 )8S0RT( AS*87)8(83-RIR)/84 )8XKU • 


ND 11*16 
ND 17 
ND 10 

ND 20 
ND 30 
ND 40 
ND SO 




sat 

C<6>«<a.*C3.14 ) »2*|7 ) 1 ( 11 -13 )/|5 >3XKU 

a 

CN 

18ft« 

S9i 

C<?>-<a.a<3.14)aSQRT(A4*a7 >*<32-11 )/|fi)«X<CU 

6 

CN 

1316 

66t 

C<a>0 3.14*<A«-A*-37 J-'IS >*XKU 

6 

CN 

1611 

611 

C<9>-<3. 14*< <*6»»7 )/3S )*XXU 

• 

CN 

1614 

6ai 

C < t • >• 3 . 1 4* < RIR-SQftT < 3 . *XN ) * 




63i 

*< <2.*«0)/4. > ttXKW 

6 

CN 

1666 

64i 

cot >03.14*<33-RIR> J*XKU 

6 

CN 

1436 

66t 

COe)*<3.14«(»l-t3)>aXKtf 

6 

CN 

1646 

661 

CO 3 >0 3. 144(82-31 ) >*XKU 

• 

CN 

1666 

671 

004 >•< 3. 14* <A7-A8-8?>/< .54(82-81 >> >*XKW 

6 

CN 

1317 

6a t 

C 05)0 3. 14«< A7-A8-87 )/< .S*< 82-81 >>>*XKW 

6 

CN 

1617 

69 1 

c oa >• c . s* < aa-at > > )«xkr 

• 

CN 

6617 

761 

co?>oa«o83'2. >>*xkr 

• 

CN 

6636 

711 

cuDxaa/oai-ffima.snaxK* 

6 

CN 

3646 

72i 

009 >0 880 (32-13 )* . 5 > )*XKR 

6 

CN 

4666 

731 

0(29)0 .714C-9*. 1*3. i4i(<4. aa( a. /3, >>* 




741 

*r!iRa*2>04.*(Ri**2>>> 

6 

CN 

1769 

751 

RETURN 



76> 

END 





EOF 176 SCAN* 16 


CO 

-1 
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WASTE, SPACE OPTION 

By Vrrv A* IVwpiriH 

Tin < mi* a* mat ion u\ ?: i< rvu<*n l .yh U'i n r\.\u%v& b > r iwhnwA 
Rtniru am inU^ ^ ;iu*h , nr;vrMr,c in i pu\. . t A IX X rauirar rwcryv M 

aetiwiu'^ Bit ill- MMM §\Mb X\ ij * ;/si C Stvwnn Cia^sitVation 

( XI HVf* 1 blS report, W- l\ % N ar 1 H/i T* ON. b. r W'liiVU X< Ik' NWi'I 



